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1.  Problems  with  boundary  layers  at  one  endpoint 

Many  physical  problems  can  be  studied  as  singularly  perturbed  two-point  vector 
boundary  value  problems  of  the  form 


where  c is  a small  positive  parameter  (cf.,  e.g.,  Amundson  (1974),  Sethna  and 
Balachandra  (1976),  and  Cohen  (1977)).  Scalar  problems  of  this  form  are  analyzed 
quite  thoroughly  in  the  forthcoming  memoir,  Howes  (1978).  An  enlightening  case 
history  of  such  analyses  was  given  by  Erdelyi  (1975),  and  Important  early  work 
Includes  that  of  Coddlngton  and  Levinson  (1952)  and  Wasow  (1956). 

For  simplicity,  let  us  assume  that  f and  g are  infinitely  differentiable 
in  y and  t and  that  they  possess  asymptotic  power  series  expansions  in  c as 
e -►  0.  We'll  first  consider  the  vector  problem  under  the  assumption  that  the 
reducfc.d  problem 


is  stable  throughout  0 ^ t jc  1 in  the  sense  that  u_  exists  and 


f is  a strictly  stable  matrix  having  eigenvalues  with  negative 
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' real  pares).  We  first  realize  that  cannot  generally  represent  the  solution 

to  (1)  near  t ••  0 because  we  cannot  expect  to  have  “ y(0).  Instead,  we 

must  expect  boundary  layer  behavior  to  occur  near  t • 0,  providing  the  required 
nonunlform  convergence  from  y(0)  to  as  £ -►  0.  For  a (very)  small 

"boundary  layer  Jump"  lly(O)  - Uj^(0)ll  or  for  a constant  f(y,0,0),  no  extra 
hypotheses  are  needed.  More  generally,  however,  we  must  require  an  additional 
"boundary  layer  stability"  assumption,  namely  that  the  inner  product 


(4) 


f(Uj^(0)+s,  0,  0)ds  > 0 


I 

remains  positive  for  c + Uj^(O)  along  all  paths  connecting  '^jjCO)  and  y(0) 
with  0 < ^ ly(0)  - Uj^(0)l . (Here,  T represents  the  transpose  and  Hz®  = 

\/z^z.)  We  note  that  if  f(z,0,0)  is  the  gradient  7F(z  - Uj^(O)),  (4)  is 

equivalent  to  the  condition  that 


C^(F(0  - F(0))  > 0 
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since  the  Integral  is  then  path- independent.  Indeed,  (4)  directly  generalizes  the 
(minim.ii  < hypotheses  used  by  Howes  for  the  scalar  problem  and  it  is  weaker  than 
the  common  assumption  that  f(y,0,0)  > 0 for  all  y. 

Pictorially,  the  boundary  layer  stability  assumption  must  hold  within  the 
circle  shown 


The  results  of  Howes  and  others  suggest  that  under  such  hypotheses,  (1)  will 
have  a solution  y(t,e)  of  the  form 


(5) 


y(t,e)  » U(t,e)  + n(T,e) 


where  the  outer  solution  U has  an  asymptotic  expansion 
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U(t,c)  -v  E U.(t)e-* 
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providing  the  asymptotic  solution  for  t > 0,  while  the  boundary  layer  correction 
n has  an  expansion 


(7) 


n(T,e)  -V  I n (T)e^ 
j-0  J 


whose  terms  all  tend  to  zero  as  the  stretched  variable 


(8) 


T » t/e 


tends  to  infinity.  We  would  expect  this  solution  to  be  unique.  Under  weaker 
smoothness  assumptions  on  f and  g,  we'd  have  to  limit  the  expansions  to  finite 
order  approximations.  For  the  scalar  problem,  Howes  doesn't  actually  obtain 
higher  order  terms  or  complete  boundary  layer  behavior,  but  they  can  easily  be 
generated.  Applying  his  results  to  the  boundary  value  problem  for  the  remainder 
terms,  however,  shows  the  asymptotic  validity  of  the  expansions  so  obtained. 

The  outer  expansion  (6)  must  provide  the  asjnnptotic  solution  to  (1)  for 
t > 0,  since  n is  then  asymptotically  negligible.  Thus,  the  terms  Uj  can  be 
successively  obtained  by  equating  coefficients  in  the  terminal  value  problem 


(9) 


f (U(t,e),t,e)U' (t,e)  + g(U(t, e) , t, e)  = -eU"(t,e),  U(l,e)  = y(l). 


Evaluating  at  e ■ 0,  then,  shows  that  must  satisfy  the  reduced  problem 

(10)  f(UQ(t),t,0)U^(t)  + g(UQ(t),t,0)  - 0,  Uq(1)  - y(l) 

(which  has  a unique  solution  u„(t)  under  (2)  and  (3)).  Succeeding  terms  U,, 


j > 0,  will  satisfy  linear  problems  of  the  form 


j' 


(11) 


f(UQ(t),t,0)Uj(t)  + fy(UQ(t),t,0)Uj(t)U^(t) 


+ gy(Ug(t),t,0)Uj(t)  - hj_j(t),  Uj(l)  - 0 

where  known  in  terms  of  t,  UQ(t),  ...,  Uj_j^(t).  The  stability  assump- 

tion (3)  implies  that  (11)  is  a nonsingular  initial  value  problem,  so  it  also  has 
a unique  solution  throughout  0 ^ t ^ 1.  Thus,  there  is  no  difficulty  in  generat- 
ing the  outer  expansion  U(t,e)  with  U(t,0)  ■ Uj^(t). 

The  boundary  layer  correction  II  must  necessarily  be  a decaying  solution  of 
the  nonlinear  initial  value  problem 
f 

2 

(12)  + f(U(cT,e)  + n(T,e),  ex,  ^ " -e[  (f  (U(ET,e)  + II(T,e),  ex,  c)  - 

dx^ 


« X 


Ifey^r'v-. . 
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- f (U(eT,e),eT,e))  ^ (er.e)  + g(U(eT,e)  + II(T,e),  ex,  e) 


(12) 


- g(U(eT,e)  ,eT,e)  ] , t >_  0 


n(0,e)  - y(0)  - U(0,e). 


Thus,  the  leading  term  11^  must  satisfy  the  nonlinear  problem 


d^n 


dn. 


(13) 


dr 


^ + f(U|,(0)  + HijCt),  0,  0)  - 0,  - y(o)  - 


“o“» 


while  later  terms  must  satisfy  linear  problems 


d^n. 


dn 


(14) 


dr 


^ + f(UQ(0)  + nQ(T),  0,  0) 


dn. 


+ fy(UQ(0)  + njj(T),  0.  o)nj(T)  ^ n^(o)  - -Uj(o) 


whei%  k.  , is  a linear  combination  of  preceding  terms  n and  their  derivatives 
j-1  *• 


dn  /dx,  i < j,  with  coefficients  that  are  functions  of  x and  ^^(x).  The 


decaying  solution  of  (13)  must  satisfy 


dn 

dx 


2+  f 1% 

* OP 


dn  • 

(0)  + n^,  0,  0)  ^ dx  - 0 


and,  thereby,  the  initial  value  problem 


dn. 


(15) 


dx 


no(x) 


f(UQ(0)  + w,  0,  0)dw,  X i 0,  ^^(0)  - y(0)  - u^(0) . 


Multiplying  by  the  boundary  layer  stability  condition  (4)  implies  that 


(16) 


n„(T) 


f(UQ(0)  + 2,  0,  0)dz  < 0 


for  nonzero  values  of  ^^(x)  satisfying  •nQ(x)l  _<  •y(0)  - Uj^(0)l  - lnQ(0)l. 

Thus,  our  boundary  layer  stability  implies  that  lnQ(x)l  will  decrease  monotonic- 
ally  as  X Increases  until  we  reach  the  rest  point  ^q(‘^)  “ ® (15)  at  x ■ 


Ultimately,  ^^^(x)  will  become  so  small  that  (3)  (for  t - 0)  implies  that  the 


eigenvalues  of  f(UQ(0)  + nQ(x),  0,  0)  will  thereafter  have  real  parts  greater 


than  some  k > 0 and  (15)  then  implies  that 


nQ(x)  - 0(e"''^), 


1 
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(17) 
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l.e.,  IIq  is  exponentially  decaying  as  x -*  «.  Although  we  can  seldom  explicitly 
Integrate  the  nonlinear  system  (15),  we  can  approximate  its  solution  arbitrarily 
closely  by  using  a successive  approximations  procedure  on  (15)  (cf.  Erdelyi 
(1964)).  Knowing  n^,  we  next  Integrate  (14)  for  j - 1 and  then  proceed  term- 
wise.  Rearranging  (14)  and  integrating,  we  obtain 


+ f(UQ(0)  + ng(T),  0,  0)nj 


+ )Ij(t) 


0 


where 


fT 


.j(T)  -= 


dn. 


{fy(Ug(0)  + nQ(r),  0,  0)[n^(r)  ~ (r) 


jf-  (r)nj(r)]  + kj_^(r)}dr 


is  known  whenever  Ilj  and  dn^/dx  commute.  Thus,  11^  satisfies  the  integral 
equation 

(18) 


n,(x)  P(x)u.(o)  - 

j i } 


P(T)P'^(r)Jlj(r)dr 


where  P(x)  is  the  exponentially  decaying  fundamental  matrix  for  the  linear  system 


+ f(u„(o)  + n,  0,  o)n  - o,  x > o,  n(o)  - i. 

dx  U 


In  general,  (18)  must  also  be  solved  via  successive  approximtlons,  though  it 
directly  provides  the  solution  of  (14)  when  the  commutator  ^IlQ/dx)  = 0. 

We  note  that  the  boundary  layer  jump  111^(0)1  •»  lly(O)  - Uj^(0)l  is  limited  by 
the  minimum  value  of  II  > 0 such  that  the  inner  product 


(19) 


f (Uj^(O)  + 2,  0,  0)dz 


0. 


The  jump  can,  in  practice,  be  quite  large.  It  Involves  no  restriction,  for  exam- 
ple, if  f(z,0,0)  > 0 for  all  z since  then  (19)  cannot  ever  hold  for  a c 0. 
It  gives  precise  limits  to  the  jumps  for  certain  scalar  problems  (Howes  (1977) 
reconsiders  an  example  of  O'Malley  (1974)). 

We  could  also  consider  the  reduced  problem 

(20)  f(Uj^,t,0)u^  + g(Uj^,t,0)  - 0, 

Then  the  stability  condition  (3)  and  the  boundary  layer  stability  condition  (4) 
would  be  replaced  by 


(22) 


f(Uj^(l)  + z,  1,  0)dz  < 0 


for  all  0 + u,  (1)  on  paths  between  u (1)  and  y(l)  satisfying  0 < 101  ^ 
L L 


ly(l)  - u (1)1 . Nonuniform  convergence  of  the  solution  to  (1)  would  then  take 


place  near  t ■ 1,  depending  on  the  stretched  variable 

« 

0 » (1  - t)/e, 


and  the  limiting  solution  on  0£  t < 1 would  be  f were  nonsingular 

with  eigenvalues  having  both  positive  and  negative  real  parts  along  an  appropriate 
solution  of  the  reduced  system,  we  must  expect  boundary  layer  beliavior  near  each 
endpoint  (cf.  Harris  (1973)  and  Ferguson  (1975)  for  discussions  of  problems  where 


fy(y.t,0)  H 0). 


2.  Problems  with  boundary  layers  at  both  endpoints 

Let  us  now  consider  the  "twin"  boundary  layer  problem 
/ 


(23) 


ey"  + R(y.t,e)  - 0,  0 £ t £ 1 

y(0),  y(l)  prescribed 


under  the  assumption  that  g is  infinitely  differentiable  in  the  region  V of 
interest  and  that  the  reduced  system 


(24) 


g(u,t,0)  - 0 


has  a smooth  solution  UQ(t)  throughout  0 ^ ^ 1 which  satisfies  the  stability 

assumption 


(25) 


gy(UQ(t),t.O)  < 0 


0 < t < 1. 


there,  i.e.,  g^  is  a stable  matrix  when  evaluated  along  (UQ(t),t,0), 

Motivation  for  this  assumption  is  obvious  if  one  considers  the  linear  scalar  prob- 


lems with  cy" ± y • 0,  while  generalized  stability  assumptions  are  sometimes 
appropriate  and  necessary  (cf.,  e.g.,  Howes  (1978)  or  consider  the  scalar  problem 
with  g • y^**^^).  With  (25),  one  can  hope  that  a solution  to  (23)  exists  which 


converges  to  UQ(t)  within  (0,1).  Since  we  won't  generally  have  either 


J 


Uq(0)  ■ y(0)  or  Uq(1)  ■ y(l)f  we  must  expect  "twin"  endpoint  boundary  layers 
(i.e.,  regions  of  nonunlform  convergence  of  thickness  0(/e)  near  both  t ■ 0 
and  t ■ 1).  Our  previous  experience  (cf.  Fife  (1973,  1976),  Yarmlsh  (1975), 
O'Malley  (1976),  and  Howes  (1978))  suggests  that  we  must  add  "boundary  layer 
stability"  assumptions.  These  generally  limit  the  size  of  the  boundary  layer 
jumps  ly(0)  - Uq(0)I  and  •y(l)  - Uq(1)I . They'll  certainly  be  guaranteed  If 
remains  stable  throughout  the  boundary  layer  regions  (cf.  Kelley  (1978)). 
Indeed,  for  small  boundary  layer  jumps,  the  stability  assumption  (25)  is  suffi- 
cient. Under  appropriate  assumptions,  then,  we  can  expect  to  obtain  an  asymptotic 
solution  to  (23)  in  the  form  , 

(26)  y(t,e)  ■ U(t,e)  + ^(D,>/e)  + w(a,/r) 

where  U,  v,  and  w all  have  power  series  expansions  In  their  second  variables 
and  the  terms  of  the  left  boundary  layer  correction  v tend  to  zero  as  the 
stretched  variable 

(27)  p - t//l 

tends  to  Infinity  while  the  right  boundary  layer  correction  w •>  0 as 

(28)  a - (1  - t)l/l 

becomes  unbounded. 

The  outer  expansion 

(29) 

should  therefore  satisfy 

(30)  eU"  + g(U,t,e)  - 0,  0 < t < 1 

as  a power  series  in  e and  converge  to  the  solution  ^^^(1)  of  the  reduced  sys- 
tem (24)  as  c ->■  0.  Higher  order  terms  In  (29)  must  satisfy  linear  systems  of  the 
form 

where  is  known  termwise  (e.g.,  ■ “U^) . The  stability  condition  (25) 

Implies  that  the  systems  (31)  are  all  nonsingular.  Therefore  successive  coeffi- 
cients are  simply  and  uniquely  obtained  termwise.  (Different  roots  Uq  of  (24) 


U(t,e)  E U (t)c^ 
j-0  3 


r 

ij 


T 


8 


would,  of  course,  result  In  different  sequences  of  perturbation  terms  , J > 0, 
under  appropriate  stability  assumptions.) 

According  to  the  Borel-Ritt  Theorem  (cf,  Wasow  (1965)),  there  is  a (non- 
unique) function  U(t,e),  holomorphic  in  e,  having  the  outer  expansion  (29). 


If  we  set 


y(t,e)  - U(t,c)  + z(t,c). 


we  convert  the  problem  (23)  into  tJje  two-point  problem 


ez"  - h(z,t,e),  d 1 t <_  1,  z(0,e)  - y(0)  - U(0,c). 


h(z,t,c)  •>  -eU"(t,e)  - g(U(t,E)  + z,  t,  e) 


satisfies 


(34)  h(0,t,e)  ••  0(e  ) for  every  integer  N >1  0 

- - M 

since  eU"  + g(U,t,c)  ■ 0(c  ).  In  particular,  the  reduced  system 

» 

h(z, t,0)  « 0 

corresponding  to  the  transformed  problem  (33)  has  the  (not  necessarily  unique) 
trivial  solution  and  the  outer  expansion  for  (33)  is  also  trivial.  Henceforth, 
then,  we  shall  deal  with  (33)  and,  corresponding  tc  (26),  we  shall  seek  an  asymp- 
totic solution  of  the  form 

(35)  z(t,E)  - v(p,/e)  + w(o,</e) 

providing  the  needed  boundary  layer  decay  to  zero  within  0 < t < 1.  Our  smooth- 
ness assumptions  will  be  required  in  a domain 

0^  g - { (z,t,£);  0 < Iz  - UQ(t)l  1 dg(t),  0 5 t 1,  0 <.  E £ E^} 

where  e^^  is  a small  positive  number  and,  for  any  6 > 0,  we  define 


- 


•z(0)  - Uq(0)*  +6,  0 < t < 6 


6, 


6 < t < 1 - 6 
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lz(l)  - Uq(1)I  +5,  l-6<t£l. 


We  shall  determine  the  asymptotic  behavior  of  z by  first  determining  that 
of  Izl  - J z^z.  Here  Izl  satisfies  the  scalar  problem 

(36)  elzl"  ■ [h^(z,t,e)z  + e(lz'l^  - (I  zl  ^)  ' ]/•  zl  , 

0 ^ t ^ 1,  where  IzCOte)!  and  llz(l,e)l  are  prescribed. 

This  follows  via  simple  calculations,'  namely 

^ Izl^  - 2«z«lzl’  » 2(z')^z 
at 


Izl^  - 2lz»«2r'  + 2(lzli  ')^ 

dt”^ 

- 2(z")^z  + 2llz'l^ 

imply  the  differential  equation  for  Izl.  Further, 

(37)  Iz’ll^  ^ (Izl  ’)^ 

T 2 2 2 2 

since  the  Cauchy-Schwarz  Inequality  ((s')  z)  ^ Iz'l  Iz#  implies  that  Hz’!  ^ 
T 2 2 

((z')  z/lzl)  ••  (Izl')  . Thus,  with  a loss  whenever  z and  z'  are  not 


collinear. 


elzl"  _>  h^(z,t,e)z/l  zl , 0 £ t £ 1. 


(Through  the  inequality  (37),  then,  we  eliminate  the  first  derivative  term  from 

(38).  We  note  that  (38)  is  an  equality  for  scalar  problems.) 

We’ll  now  ask  that  for  all  (z,c,e)  in  P ,,  there  exists  a smooth  scalar 

e^,6 

function 


♦(n,t,e) 


such  that 


h (z,  t,e)z  3^  ♦ (I  z* , t,e)l  zl 


I 


where 


whenever  0 < n < 


whenever 


Existence  of  such  a function  ^ will  constitute  our  stability  hypotheses.  Spe- 
cifically, (0,t,0)  > 0 implies  the  stability  of  the  trivial  solution  of  the 

reduced  system  within  (0,1)  while  (Al)  implies  boundary  layer  stability  at  l>oth 
endpoints.  Hypotheses  (39)- (40)  imply  that 


where  m(t,c)  satisfies  the  scalar  two-point  problem 


The  bounds  (42)  follow  from  the  elemei\tary  theory  of  differential  inequalities 

since  zero  is  a lower  solution  for  I zl  and  m is  an  upper  solution  (cf.  Nagumo 

(1937),  Dorr,  Parter,  and  Shampine  (1973),  and  Howes  (197b)).  Further,  ^(0,t,0) 

■ 0 and  (0,t,0)  > 0 imply  that  the  zero  solution  of  the  reduced  problem 
3n 

<^(m,t,0)  0 corresponding  to  (43)  is  stable  and,  according  to  Howes  (1978),  (41) 

is  the  appropriate  hypothesis  for  the  needed  boundary  layer  stability  of  this 
solution.  Indeed,  the  solution  of  (43)  satisfies 


where  r is  the  decaying  solution  of  the  boundary  layer  problem 


is  the  decaying  solution  of 


while  s 


The  solutions  to  (A5)  and  (A6)  are  easily  shovm  to  exist  and  be  unique.  Multiply 
ing  (45)  by  dr^/dp,  for  example,  and  integrating  from  p to  infinity  implies 
that 

^(s,0,0)ds  > 0 


(by  (41)).  Thus, 


satisfies  the  initial  value  problem 


(47) 


Hence,  will  decrease  monotonically  to  zero  as  p increases,  reaching  the 

rest  point  r.  ■ 0 at  p ■ “>.  Since  iji(s,0,0)  (0,0, 0)s  for  s small, 

g.  0 3n 

^ (0,0,0)  > 0 implies  that  the  decay  of  r.  to  zero  is  exponential  as  p 
3n  u 

(When  *'q(0)  “ 0*  have  ® since  there  is  no  need  for  a boundary 

layer  correction.)  Continuing  by  solving  linear  problems,  we  could  obtain  an 
asymptotic  solution  of  (43)  in  the  form 


(48) 


m(t,c)  - r(p,/c)  + s(o,/e). 


In  terms  of  the  original  problem  (23),  our  stability  hypothesis  (39)  becomes  the 
inequality 


g^(UQ(0)  + z,  t,  e)z  £ (}i(l  zl  ,t,e)i  zl 


where  ^ satisfies  (40)  and  (41).  The  expansion  (44)  corresponds  to  the  expected 
expansion  (26)  for  an  asymptotic  solution  for  the  vector  problem  (23). 

Now,  we  return  to  the  vector  boundary  value  problem  (33)  and  its  asymptotic 
solution  in  the  form  (35).  Near  t - 0,  w and  its  derivatives  should  be  asymp^ 
totlcally  negligible  (o  being  infinite),  so  (33)  and  (35)  imply  chat  the  initial 
boundary  layer  correction  v should  be  a decaying  solution  of  the  nonlinear 
Initial  value  problem 

(49)  V ■ h(v,/ep,e),  p > 0,  v(0,e)  ■ z(0,e). 

PP  — 

Thus,  it  is  natural  to  seek  an  expansion 

(50)  v(p,r^)  'v-  I V (p)e^^^ 

J-0  J 

by  substitution  into  (49).  The  leading  term  v^  must  then  satisfy  the  nonlinear 


problem 


12 


(51)  — ^ - hCv^.O.O),  P > 0,  Vq(0)  - y(0)  - Uq(0).  0 as  p 

dp 


Later  terms  , j i 1,  must  satisfy  the  linear  problems 


- h (v  ,0,0)v  + d ,(p),  p ^ 0 

2 z 0 j J-1  - 


(52)  \ Vj(0)  - 0,  j odd;  Vj  (0)  - ^ 

V -^0  as  p -*■  “> 

where  d^  ^ will  be  determined  successively  as  an  exponentially  decaying  vector, 
.'liiice  (5i)  and  our  hypothesis  (39)  imply  that 


“Vq»pp  i <^(Ivq«,0,0),  p > 0,  »Vq(0)I  - rQ(0), 


we  are  guaranteed  a decaying  solution  such  that 


0 £ Jvq(p)1  _<  TqCp).  P 1 0 


(and  ''q(p)  =0  ” y(0)).  No  explicit  solution  v^  can  be  provided, 

though  an  approximate  solution  can  be  obtained  as  usual.  Introducing  tlie  matrix 

a = h (0,0,0)  > 0 
z 

(whose  eigenvalues  have  strictly  positive  real  parts  by  our  stability  assumption 
(25)),  variation  of  parameters  can  be  used  to  express  the  solution  of  (52)  in  the 


V (p)  - e V (0)  - [ e (s)ds 

J J J Q J _ 

1 / /“\“1  -»^(p-r)„  , ,,  , f *^(p-s) 

- — (/a)  e F,(r)dr  + e F (s)ds 

2 [Jo  J Jp  J 


where  F (p)  ■ [h  (v  (p),0,0)  - a]v  (p)  + d ,(p).  This  provides  the  exact  solu- 
J z 0 J j“l 

tlon  to  (52)  whenever  h(vQ,0,0)  is  linear.  Otherwise,  the  linear  integral 
equation  (54)  must  also  be  solved  by  successive  approximations.  In  analogous 
fashion,  we  could  generate  the  terms  of  the  terminal  boundary  layer  correction 
w(o,/e)  of  (35).  Thus,  we've  formally  obtained  (35),  which  we  expect  is  a locally 


unique  asymptotic  solution. 


We  note  that  the  assumptions  on  automatically  hold  if  gy(y.t,0)  or 
-h^(z,t,0)  are  everywhere  stable.  Thus,  if  we  take 

h (v,0,0)  - Yl  > 0 
z 

(i.e.,  positive  definite)  for  some  real  y > 0 and  all  v satisfying  0 <.  Ivl  £ 
llz(0,0)ll,  the  mean  value  theorem  implies  that 

h^(z,0,0)z  ••  z^h  (z,0,0)z  > 
z 

I 

for  some  "intermediate"  point  z.  Thus,  taking  il>(n,0,0)  - yn,  both  (0,0,0) 

> 0 and  /q  4i(s,0,0)ds  > 0 for  0 < n f.  llz(0,0)ll  hold. 

We  could  also  extend  our  discussion  to  systems  of  the  form 


ex"  - F(x,x' ,t,e) 


3F 


with  3F/3x'  small.  Thus,  Kelley  (1978)  considered  problems  where  ~ 

_L  IL.  IL.  I > 0,  just  as  Erdfilyi  (1968)  considered  scalar  problems  somewhat 
2c  3x  \ 3x  / 

more  nonlinear  than  semllinear. 


3.  Examples 

a.  A problem  with  an  initial  boundary  layer 
Let  us  consider  the  vector  equation 

ey"  + f(y,t,e)y'  + g(y,t,e)  - 0,  0 £ t £ 1 


where 


1 

1 h 1 1 

f * M 

, f - 

, and  g “ - 

^2  i 

>'2i 

i ^2  1 1) 

In  order  to  have  a limiting  solution  u„  of  the  two-point  problem  which  satisfies 

K 

the  reduced  problem 

f(uj^,t,o)u^  + g(u^,t,o)  ■ 0, 

we  must  require  u^^  to  be  stable  in  0 £ t £ 1,  i.e., 

-f(Uj^(t),t,0)  < 0 


must  be  a stable  matrix,  and  we  must  also  require  boundary  layer  stability  at 


J 


(We  note  that  the  minimum  for  “ 0.  t - is  2D.)  This  calculus  problem, 

then,  determines  an  upper  bound  for  ly(0)  - u„(0)l. 

3 ^ 

For  C - D - 2,  i.e.,  y(l)  - (p,  we'd  obtain  the  minimum  value  3.390 

for  d(t)  corresponding  to  t^^^  - -0.291.  Thus,  we're  guaranteed  that  the 

limiting  solution  of  our  two-point  problem  is  provided  by  if  y(0)  lies 

2 

in  the  circle  of  radius  3.390  about  C^).  This  is  presumably  a conservative 
estimate  for  the  "domain  of  attraction"  of  the  reduced  solution  We  expect 

that  boundary  layer  stability  need  only  hold  for  C + actual  trajec- 

tory joining  y(0)  and  Uj^(O).  Finally,  we  observe  that  this  example  is  quite 
analogous  to  the  simplest  cases  occurring  in  the  analysis  of  solutions  of  the 
scalar  problem  ey"  + yy'  ~ y ■ 0 (cf.  Cole  (1968),  Howes  (1978),  and  elsewhere). 

b.  A problem  with  twin  boundary  layers  at  the  endpoints 
Consider  the  vector  problem 


where 


2 » 


the  Jacobian  matrix 


ez" 

■ h(z,t,e). 

0 <.  t <.  1 

' 'l  * '2  • \ 

^ "2  i 

and  h ■ 

. -'l  *^2-4 1 

e solution  of  the  reduced  problem 

1 1 1 \ 

h^(0,t,0)  - 

-1  1 j 

0 since 


has  the  unstable  eigenvalues  1+1.  Boundary  layer  stability  involves  the  deter- 
mination of  a scalar  function  ^ such  that 


h^(2,t,e)z  ^ ^ (I zl , t ,e)l zl . 


Here 


h^(z,t,e)z  - (z^  z^)  - (zj  + z^)  >1  Izi^d 


• zl^). 


^ ^ ^ / 2 ^ 2.2 

Since  Zj^  + z^  <.  (z^  + z^)  , so  we  can  take 


4>(n,t,c)  • n(l  - n ). 


Clearly,  ^(0,t,c)  <.  0,  ♦(0,t,0)  - 0,  ♦^(0,t,0)  > 0 and 


mrnmmmm 
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(;>(s,i,0)ds  •■  -r  n (1  - n /2)  >0  for  0 < n < >'2,  i - 0 or  1. 

Jo 

Our  preceding  results,  then,  guarantee  the  existence  of  an  asymptotic  solution  to 
the  two-point  problem  which  converges  to  the  limiting  solution  0^-0  within 
(0,1)  provided  the  boundary  values  satisfy 

llz(0,0)ll  < /2  and  llz(l,0)ll  < /I, 


Indeed,  we  then  have 


0 £ II  z(t,e)ll  _<  m(t,e) 

where  m satisfies  the  scalar  problem 

em"  - ((i(m,t,e),  0 £ t £ 1,  m(i,e)  - llz(l,e)ll  < >^,  1-0  and  1. 

The  asymptotic  behavior  of  m follows  from  the  scalar  results  of  Howes  (1978) 
and  others. 

c.  A problem  with  internal  ir.insitlon  layers 
We  now  consider  the  very  special  problem 


1 


ey"  + f(y,t,e)y'  + g(y,t,e)  - 0,  0 £ t £ 1 

where 


f h } 

1 

^2^yi’y2’^’^J  ' 

y ■ 

. f(y.t,e)  « 

{ ^2  1 

\ 0 

yz 

and  g ■ 

\ -y2 


This  system  decouples  into  the  two  nonlinear  scalar  equations 


ty'^  + y2y2  ■ y2 


0 


ey'^  fj^(yj.t,e)y|  + 


0. 


and 
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If  y2(l)  > 72(0)  + 1 and  -y2(l)  - 1 < y2(0)  ^ ^ follows  from 

Howes  (1978)  that  the  limiting  solution  for  y^  will  satisfy  the  reduced  problem 

Uj^(u^  - 1)  - 0,  Uj^(O)  - y2(0)  on  0 < t < t*  - ^ ^1  - y2(l)  - y2(0)) 

and  the  reduced  problem 


1 -t 


1 


Uj^(u^  - 1)  - 0,  Uj^(l)  - on  t*  < t < 1, 


i.e. , 


u- 


u (t)  “ t + y„(0),  0 t < t* 


Uj^(t)  - t + y^Cl)  - 1,  t*  < t £ 1. 


Thus,  the  limiting  solution  is  generally  discontinuous  at  t*  and  its  derivative 
(which  is  asymptotically  one  elsewhere)  becomes  unbounded  there.  Indeed,  y^ 
increases  monotonically  near  t*  from  u^(t*)  to  Uj^(t*)  *■  -Uj^(t*).  For  other 
relations  between  the  boundary  values  y2(0)  and  y2(l),  other  limiting  possi- 
bilities occur  (cf.,  e.g.,  Howes). 

One  must  generally  expect  the  transition  layer  at  t*  in  y2  to  generate  a 
corresponding  discontinuity  there  in  y^^.  To  simplify  our  discussion,  however, 
let's  assume  that  f2(yj^,y2.t.0)  * 0 ^nd  attempt  to  apply  Howes'  scalar  theory 
to  the  equation  for  y^^.  Thus,  consider  the  reduced  problems 

' gj(Vj^,u,t,0)  - 0,  0 <.  t <.  1,  “ yj^^o) 


and 


fl(vR,t,0)v^  + 8j^(vR,u,t,0)  - 0,  0 £ t £ 1,  “ y2^(l)- 

The  limiting  solution  for  y^^  will  be  provided  by  ^^(t)  if  the  stability  con- 
dition 


fl(vR(t),t,0)  > 0 


holds  throughout  ^ ^ f ^ f and  the  boundary  layer  stability  assumption 

/r“> 

(Vr(0)  - yj^(O))  fj^(s,0,0)ds  > 0 

* n 


for  n between  Vo(0)  and  (including)  y, (0).  Similar  conditions  would  imply 
K X 

that  the  limiting  solution  is  Vj^(t)  on  0 £ t < 1 with  boundary  layer  behavior 


Note  that  has  a jump  at  t and  has  a jump  at  t*,  corresponding  to  a 
Haber-Levinson  crossing  (cf.  Howes  (1978)),  Much  more  complicated  possibilities 
remain  to  be  studied. 
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